Mid and Far Infrared Surveys with ISO 



D. Elbaz \ 2 

^ CEA Saclay/DSM/DAPNIA/Service d'Astrophysique, Orme des Merisiers, F-91191 



Q\ I Gif-sur-Yvette Cedex, France 

^ ■ ^ University of California, Santa Cruz (U CSC), Physics Department, Santa Cruz, CA 95064, 

> USA 
O ■ 

^: 

■ Abstract 

I We present the results of the extragalactic surveys performed at 15 and 175/imusmg 

f-*) ' the ISOCAM and ISOPHOT instruments on-board the Infrared Space Observatory (ISO). 

if^ i The number counts at both wavelengths present a strong evolution. The preliminary study 

' of the individual galaxies responsible for this evolution suggest that more star formation 

. was hidden by dust at redshifts around z = 1—2 and imply that the population of luminous 

Q>^ I IR galaxies detected by IRAS in the local universe must have played a major role in the 

0^ ■ formation and evolution of galaxies and cannot be considered anymore as extreme and 

^ I rare cases. 

O 



■ 

X 



1 Introduction 

Several critical observational advances have taken place in the past few years to elucidate the 
formation and evolution of galaxies and the rates at which stars formed at different cosmological 



^ epochs. With the advent of magnitude limited samples of spectroscopic redshifts, it became 
. possible to discuss general properties of galaxies through their comoving luminosity density or 
star formation rate (SFR) as a function of redshift (Lilly et al 1996, Madau et al 1996). Since 
then, we became familiar with the study of a new astronomical character, the '1 Mpc^' box 
of the universe. At the same time, the plot of the cosmic background was enriched by COBE 
(Puget et al 1996, Fixsen et al 1998, Lagache et al 1999, Hauser et al 1999, Dwek et al 1998) and 
the HST (Pozzetti et al 1998), on the far infrared to sub-millimeter part and on the optical-UV 
part respectively. These two global views of the past history of galaxy formation and evolution 
gave a new life to the debate that begun with the discovery by IRAS of Luminous InfraRed 
Galaxies (LIRGs, see the review by Sanders & Mirabel 1996). Galaxies brighter than 10^^ Lq 
radiate most of their light in the infrared (IR) range and exhibit the strongest star formation 
rates ever measured. No other observations, else than direct infrared observations, had prepared 
such a discovery. However, the infrared radiation of nearby galaxies {z < 0.2) makes only 30 % 
of their optical light and only 6 % of it comes from LIRGs (based on IRAS results, Soifer & 
Neugebauer 1991), which consequently weight only 2 % by comparison with the optical light 
from 'normal' galaxies. Prior to the COBE, SCUBA and ISO results, hence only three years 
ago, the matter of debate could have been summarized by the following question: should we 
consider LIRGs, as anything else than extreme cases, non representative of the general behavior 
of galaxies ? 



When the previously quoted authors, after Puget et al (1996), found from the COBE FIRAS 
and afterwards DIRBE measurements, that the energetic content of the IR to sub-mm cosmic 
background was as high as that of the UV to near-IR background, the debate became more 
intense but still left opened the possibility to keep a scenario consistent with the one deduced 
from the '1 Mpc^' box of the universe story, i.e. that one could reproduce both the CIRB 
and the variation of the comoving UV luminosity density and SFR with redshift without the 
recourse of LIRGs (see Fall, Chariot & Pei 1996). 

Several observations have now ruled out this picture. The redshift dependance of the co- 
moving star formation rate of the universe needs not only the UV luminosity density to be 
known as a function of redshift but also the IR luminosity density, which is dominated by 
LIRGs at z ~ 1 and above (see below). The most directs of these observations, because they 
measure the IR light radiated by dust, come from three instruments: ISOCAM (Cesarsky et al 
1996), in the mid-IR, and ISOPHOT (Lemke et al 1996), in the far-IR, on board ISO (Kessler 
et al 1996), and the SCUBA bolometer array at JCMT, in the sub-mm (see Lilly et al, in these 
proceedings). All three instruments present number counts much above expectations from no 
evolution models based on IRAS results on the local universe, indicating that the bulk of the 
cosmic IR background (CIRB) originates from LIRGs, or even brighter galaxies. Hence, an- 
swering NO to the previous question. So that we can say that we are leaving in the middle of 
an "Infrared Revolution". The following discussion will illustrate the role that the ISOCAM 
and ISOPHOT extragalactic results are playing in it, being complementary to and consistent 
with the SCUBA results. 

2 Description of the Surveys and Number Counts 

ISOCAM is a mid-IR camera which allows to perform either wide band imagery or low resolution 
spectroscopy using Circular Variable Filters, in the range 5-18 /im. The two wide band filters 
that were used for the extragalactic surveys, LW2 & LW3, are centered at respectively 6.75 /im 
(5-8 fim) and 15 /im (12-18 fim) and were chosen because they cover the mid-IR emission from 
two different origins, the aromatic features and the hot dust continuum (see Section |^). Its 
pixel size is 6" in the most sensitive mode, for a FWHM PSF of 9" at 15 /im, but its astrometric 
accuracy can reach 2" in the micro- scanning mode of the deepest surveys. 

ISOPHOT is an imaging photo-polarimeter covering the 2.5-240 fim range. The most per- 
formant filter for extragalactic surveys is centered at 175 yum (130-219yum) and uses the four 
pixels camera, C200, with a 92" side length and a 1.9' FWHM PSF at 175 /im. 

The list of all ISOCAM extragalactic surveys performed during the ISO lifetime is given 
in the Table ^ We used 8 of them to produce the number counts in the Fig. |^ (from Elbaz 
et al 1999), for a total of 614 sources above an 80 % completeness limit. The ISOPHOT 
number counts from Dole et al (1999) superimposed in the Fig. ^ (left) use 208 sources brighter 
than 100 mJy in ~3 square degrees from the fields quoted in the Table ^. The ISOCAM 
and ISOPHOT surveys were performed in regions selected for their low Galactic foreground 
emission, both in the northern and southern hemispheres to avoid strong contamination from 
large-scale structure. 

3 Origin of the infrared emission 

The origin of the mid-IR emission of galaxies has been a subject of debate prior to IRAS, when 
evolved stars such as OH/IR stars were expected to play a major role. However, based on the 



Table 1: Table of all ISOCAM extmgalactic surveys sorted by increasing depth. The results 
used for the number counts are denoted by a -k. Col. (4) is the integration time per sky position. 
The depth of each survey is given in Col. (5), where the number between parenthesis correspond 
to the 80 % completeness limit, and the asscociated numbers of sources are given under Col. 
#7 and #15 for the 7 and 15 /im bands respectively. The coordinates are in J2000. References: 



a-in preparation, 1-Siebenmorgen et al (1996), 2-Rowan-Robinson et al (1999), 3-Clements et 
al (1999), 4-Flores et al (1999a,b), 5-Metcalfe et al (1999), 6-Aussel et al (1999), 7-Taniguchi 
et al (1997a), 8-Taniguchi et al (1997b), 9-Altieri et al (1999). 
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Tabic 2: Table of ISOPHOT-FIRBACK extragalactic surveys. References: 1-Oliver et al, in 
preparation, 2-Dole et al 1999, Puget et al 1999. 
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IRAS Point Source Catalog, source counts at 12 /xm have shown that circumstcllar dust shells 
accounted for only ~10% of the 12 yum emission in the Galaxy at | 6 |> 10° (see Review by 
Cox & Mezger 1989). It is now established that except for ellipticals (with no dust lanes), 
the mid-IR emission of galaxies does not originate primarily from individual stars but rather 
from the light of young stars (or of an AGN) absorbed by dust and re-radiated in the mid- 
IR mostly at the interface between HII regions and molecular clouds (the photo-dissociation 
regions, PDR, see Cesarsky & Sauvage 1999 and references therein). The mid-IR spectrum of 
a galaxy can be divided into three components: a set of broadband aromatic features, with 
their underlying continuum, narrow forbidden lines of ionized gas (negligible in ISOCAM wide 
band imagery), and a hot thermal continuum due to very small grains (VSGs) of dust heated 
by stellar radiation at temperatures greated than 150 K. The broadband aromatic features 
have been proposed to originate either from PAHs (for Polycyclic Aromatic Hydrocarbons, 
Lcger & Puget 1984, Allamandola et al 1989) or from coal grains (Papoular 1999). In either 
case, they require the presence of aromatic structures, but their exact nature has not yet been 



definitively established. The low resolution spectra from ISOCAM have clearly established 
that in the wavelength range 5-8.5 /xm of the 7/im filter, the spectra of spiral and starburst 
galaxies were dominated by the aromatic broadband features located at 6.2, 7.7, 8.6, 11.3 and 
12.7yum. At redshifts of the order of z~0.8, the median redshift of the faint ISOCAM sources 
(see Section |^ ), the 15 /xm band is shifted to this region of the spectrum dominated by the 
aromatic features, with a contribution of hot dust continuum. 

At larger wavelengths, the spectral energy distribution (SED) of a galaxy is progressively 
dominated by the thermal emission of larger grains heated at a lower temperature (around 
60 K), which produces the bulk of the overall IR emission originating from dust (around 60- 
100 yum). This 'cold' dust dominates the IR emission measured by ISOPHOT. At 175 /xm, the 
physical origin of the emission remains the same up to high redshifts and is well-known to be 
correlated with star formation and radio activity (Helou, Soifer & Rowan- Robinson 1985). 

The link between mid-IR and star formation (SF) is less straightforward, but only LIRGs 
can be detected at z ~ 0.8 above 0.1 mJy at 15 fim. Hence mid-IR surveys can at least pinpoint 
IR active galaxies. Moreover, a correlation was found between mid-IR and UV flux densities 
(Boselli et al 1997, 1998) as well as with (for the disks of spiral galaxies, Roussel et al, these 
proceedings). However, the complex link between mid-IR and SF is illustrated by the Fig. ^ 
(left), where we compare the SED of the ultra-LIRG Arp 220 and the nearby starburst M82, 
normalized by a factor 5 to reach the same energy density at ~8 fim. At z~0.8, both galaxies 
would have the same 15 fim flux density but differ by a factor ~20 in the far-IR. But Arp220 
is an extreme case which can easily be rejected as a candidate SED for the ISOCAM galaxies 
(see Section ^^) , and M82's SED, which is more typical of most LIRGs {Lboi ~ Ljr > 10^^ Lq) 
and ultra-LIRGs {L^oi > 10^^ Lq) although it is only 4 10^° Lq, is probably a better candidate. 
Indeed, in the HDF-N seven galaxies, which were detected both at 15 fim and in the radio 
(Richards et al 1998), show consistent SF rates if we use the SED of M82 to convert the mid-IR 
flux into a SF rate (Aussel et al 1999, in preparation). 




Figure 1: (left:) Comparison of the Arp 220 and M82 SEDs. (right:) I(Kron-Cousins)-(H+K) 
color (from Barger et al 1999) versus redshift for field galaxies (dots) and galaxies detected 
above 0.1 mJy at 15 ^m by ISOCAM (dots surrounded by circles), in the Hubble Deep Field 
North plus its Flanking Fields. Plain line: linear fit to the colors of ISOCAM galaxies. Dashed 
line: idem for all galaxies. 



4 Discussion about the 15 fim. results 



4.1 The 15 /im counts 

The first striking result of tlie 15 fim source counts is tlie consistency of the eight surveys (noted 
with a-k in the Table |1|) over the full flux range (Fig. 0). Some scatter is nevertheless apparent; 
given the small size of the fields surveyed, we attribute it to clustering effects. The two main 
features of the observed counts are a significantly super-euclidean slope (a = —3.0) from 3 
to 0.4 mJy and a fast convergence at flux densities fainter than 0.4 mJy. In particular, the 
combination of five independent surveys in the flux range 90-400 /iJy shows a turnover of the 
normalized differential counts around 400 /iJy and a decrease by a factor ~ 3 at 100 /xJy. We 
believe that this decrease, or the flattening of the integral counts, below ~400 /xJy, is real. 
It cannot be due to incompleteness, which was quantified using the Monte-Carlo simulations. 
The differential counts can be fitted by two power laws by splitting the flux axis in two regions 
around 0.4 mJy. In units of mJy~^ deg~^, we obtain, by taking into account the error bars {S 
is in mJy): 



dNjS) 
dS 



(2000 ±600) ... 0.1<5<0.4 

(1) 

(470 ±30) 5(-3-o±o-i) ... 0.4<^<4 



The total number density of sources detected by ISOCAM at 15 fim down to 100 /iJy (no 
lensing) is (2.4 ± 0.4) arcmin"^. It extends up to (5.6 ± 1.5) arcmin"^, down to 50 yuJy, 
when including the lensed field of A2390 (Altieri et al 1999). The differential counts (Fig. ^ 




Figure 2: (Left:) Integral counts, i.e. the number of galaxies, N, detected at 15 and 175 
(from Dole et al 1999, normalized in vSy) above the flux S(mJy), with 68 % confidence contours. 
K counts (Gardner et al 1993) and B counts (Metcalfe et al 1995), multiplied by the ratio 
z//z/i5 to represent the relative energy densities, are superimposed with open circles and filled 
squares, respectively. The hatched area materializes the range of possible expectations from 
models assuming no evolution at 15fim and normalized to the IRAS 12fim local luminosity 
function (LLP). (Right:) Differential Number Counts of the 15fim Galaxies, with 68% error 
bars. The counts are normalized to a Euclidean distribution (a = —2.5). 



(right)), which are normahzed to S~'^'^ (the expected differential counts in a non expanding Eu- 
chdean universe with sources that shine with constant luminosity), present a turnover around 
Si5^m=0.4: mJy, above which the slope is very steep (a = —3.0 ±0.1). No evolution predictions 
were derived assuming a pure ^-correction in a ffat universe (go = 0.5), including the effect of 
the aromatic features in the galaxy spectra. The lower curve of the hatched area, which materi- 
alizes the 'no evolution' region in the Fig. |^, is based on the Fang et al (1998) IRAS 12 fim local 
luminosity function (LLF), using the spectral template of a relatively quiescent spiral galaxy 
(M51). The upper curve is based on the Rush, Malkan & Spinoglio (1993) IRAS-12/im LLF, 
translated to 15 /im using as template the spectrum of M82, which is typical of most starburst 
galaxies in this band. More active and extincted galaxies, like Arp220, would lead to even lower 
number counts at low fluxes while flatter spectra like those of AGNs are less fiat than the one 
of M51. In the absence of a well established LLF at 15 fim, we consider these two models as 
upper and lower bounds to the actual no-evolution expectations; note that the corresponding 
slope is ~ —2. The actual number counts are well above these predictions; in the 0.3 mJy to 
0.6 mJy range, the excess is around a factor 10: clearly, strong evolution is required to explain 
this result (note the analogy with the radio source counts, Windhorst et al 1993). 

For comparison, we have superimposed on Fig. ^ the optical B band (Metcalfe et al 1995) 
and near-IR K band (Gardner et al 1993) integral counts, normalized in terms of uS^. For 
bright sources, with densities lower than 10 deg~^, these curves run parallel to an interpolation 
between the ISOCAM counts presented here and the IRAS counts; the bright K sources emit 
about ten times more energy in this band than a comparable number of bright ISOCAM sources 
at 15 fim. But the ISOCAM integral counts present a rapid change of slope around 1-2 mJy, 
and their numbers rise much faster than those of the K and B sources. The sources detected 
by ISOCAM are a sub-class of the K and B sources which harbor activity hidden by dust. 

4.2 Nature of the 15 fim galaxies 

We believe, according to the results obtained on the HDF (Aussel et al 1999a,b) and CFRS 
fields (Flores et al 1999a,b, see also these proceedings), that the sources responsible for the 
'bump' in the 15 fim counts are not the low mass faint blue galaxies which dominate optical 
counts and have a median redshift around z~0.6 (Pozzetti et al 1998). Instead, they most 
probably are bright and massive galaxies whose emission occurs essentially in the IR and could 
account for a considerable part of the star formation in the universe at z < 1.5. Indeed, from 
the full sample of 15 fim galaxies with known redshift and optical-near IR magnitudes, we found 
that these galaxies are massive (~ 10^^ ^^0) and that their emission occurs essentially in the 
IR. Their median redshift is z~0.8 in a sample of 42 galaxies brighter than 0.1 mJy in the 
HDF+FF (Aussel et al 1999b) and z~0.7 in a sample of 41 galaxies brighter than 0.35 mJy 
in the CFRS-1415 field (Flores et al 1999b). The assessment of their bolometric luminosity 
requires the assumption of a spectral energy distribution, which is largely uncertain since the 
ratio of the far-IR over mid-IR flux densities is highly variable among galaxies. However, one 
can set limits to the bolometric luminosity using the following arguments. If all galaxies had 
extreme SEDs like the one of Arp220, they would produce a contribution to the SCUBA-850 fim 
number counts (see Blain et al 1999) much above observations (see Fig. ^). An Arp220-like 
SED can also be rejected for over-producing the 140 /xm EBL with respect to the one measured 
by DIRBE on-board COBE (see Fig. ^). The light radiated in the far-IR by these galaxies 
however cannot be much lower than that emitted at 15 /im, since only galaxies luminous in 
the IR can be detected at z>0.5 with a 15/im flux density larger than 0.1 mJy (the origin of 
the emission cannot be stellar without requiring excessive masses). We therefore estimate the 
mean bolometric luminosity of these galaxies to be of the order of a few 10^^ Lq. 
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Figure 3: Integral counts at 850 [im from SCUBA (Blain et al 1999). From the 15 [im ISOCAM 
counts (dotted lines), one can extrapolate the contribution of the ISOCAM galaxies at 850 ^m 
for a given redshift distribution and spectral energy distribution (SED). Two examples are shown 
here using the SEDs of M82 and Arp 220, and assuming that all ISOCAM galaxies have the 
same redshift, z=0. 1, corresponding to their median redshift in the HDF and CFRS fields. 

The galaxies from the HDF-N plus Flanking Fields whose 15 /im flux density is greater than 
0.1 mJy (sensitivity limit of ISOCAM) harbour a Ic-(H+K) color distribution very similar to 
that of the whole sample of galaxies for which we have access to both the Ic-(H+K) colors and 
spectroscopic redshifts (see Fig. |I] (right)). Aussel et al (1999b and PhD thesis) showed that 
the sub-sample of galaxies with known spectroscopic redshift keeps the same color properties 
than the full sample of HDF+FF galaxies, hence we did not include a selection bias in the color 
distribution. 

The median colors are Ic-(H-|-K)= 2.3 for the 15 /im galaxies and Ic-(H-|-K)= 2 for all 
galaxies, hence only differ by 0.3 dex. However, the linear fit to the Ic-(H-I-K) versus redshift 
plot of the two samples of galaxies in Fig. |] (right) shows that the 'dusty' galaxies tend to redden 
slightly faster with increasing redshift, than the natural reddening of the whole population of 
field galaxies which is due to /^-correction: the color difference increases from 0.2 dex below 
z=0.7 to 0.4 dex above z=0.7. Nevertheless, this difference is not strong enough to allow one 
to separate the infrared galaxies from the whole sample like for Lyman-break galaxies (even 
accounting for the other set of optical colors existing for the HDF galaxies). The origin of this 
reddening with redshift is not clear but it is probably due to a selection of the galaxies suffering 
from more extinction, hence emitting more in the infrared, at larger redshifts. 

Hence the population of galaxies producing the 15 fim number counts excess is very distinct 
from the one which dominates the deep optical counts, known to be made of low mass galaxies 
with blue luminosities, although at similar redshifts. In other words, the star formation activ- 
ity responsible for the light emitted by the 15 fim galaxies is not the other face of the same 
star formation activity already quantified from the optical surveys, but instead it should be 
considered as a second component, which was previously missed. 
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Figure 4: Cosmic background from UV to mm. Open squares give the lower limits from ISO 
(ISOCAM-15 ixm and ISOPHOT-175 fim) and SCUBA-850 fim (the two thick crosses corre- 
spond to the recent values at 450 and 850 ^m from Blain et al 1999). The optical- UV points are 
from Pozzetti et al (1998). The COBE FIRAS (grey area) and DIRBE I40 & 24O fim (filled 
circles) data are from Lagache et al (1999). 

5 Discussion about the 175 fim results 

5.1 The 175 iim counts 

As for the 15 fim ISOCAM counts, the 175 /xm ISOPHOT counts present a strong excess with 
respect to IRAS of about one order of magnitude (Puget et al 1999, Dole et al 1999). The slope 
of the integral counts is a = —2.2, hence much above the Euclidean value of -1.5, and far from 
flattening, contrary to the ISOCAM counts which go deep enough to detect the beginning of 
the flattening, with a slope decreasing from -2.4 around 0.5 mJy to -1.2 around 0.1 mJy. It is 
therefore natural that the detected sources account only for 10 % of the CIRB measured by 
COBE (see Fig. §. 

5.2 Nature of the 175 /im galaxies 

The search for optical counterparts to the ISOPHOT-175 fim galaxies is not an easy task because 
of the large uncertainty on the position of the ISOPHOT sources and of the absence of a 
clear optical signature. Hence the use of other wavelengths with better spatial resolution 
and sensitivity to the IR activity is mandatory in this study. Because of the good radio/far- 
IR correlation, the radio domain is an obvious candidate and VLA surveys of the northern 
fields (ELAIS Nl & N2) have been carried out by Ciliegi et al (1999) in the framework of 
the ELAIS project (Rowan- Robinson et al 1999). A sample of ten 21 cm radio counterparts 
of the ISOPHOT galaxies, i.e. with Sn^^m > 100 mJy, have been selected for their high 



175yum/21 cm ratio by Scott et al (1999) for a follow-up at 450 and 850 yum with SCUBA. 
Assuming a dust temperature of = 40 K, they find a range of low to moderate redshifts 
of 2; = — 1.5 (no spectroscopic redshift are available yet), when accounting for the SCUBA 
measurements or upper limits together with the ISOPHOT ones. 

In these two fields, the ISOCAM surveys are not sensitive enough to strongly constrain 
the SED of the ISOPHOT galaxies and miss the majority of them (5*15 > 3 mJy). In the 
southern field however (called MFB, for 'Marano FIRBACK' since it is shifted with respect 
to the original Marano Field, Marano, Zamorani & Zitelli, 1988, in order to minimize the 
foreground Galactic cirrus emission), deep ISOCAM surveys have been performed over 0.5 
square degree down to 0.3 mJy (and a completeness limit of 0.4 mJy). Half of the 24 FIRBACK 
galaxies in the field have a counterpart, sometimes several, at 15 /xm above 5*15 = 0.3 mJy 
and closer than 30" (except one at 50"). One can compare their Sn^^m/ Si^^m ratio to the 
one of any template SED at any redshift and determine a possible redshift range for this SED. 
The striking result of this study is that none of the 12 galaxies present any 15 fim counterpart 
compatible with the SED of M82 at any redshift (unless they are brighter than 10^^ Lq), while 
all but one of them are consistent with an Arp220-like galaxy of similar luminosity (~ 10^^ Lq) 
and at a redshift between z=0. 1-0.4. For the remaining 12 galaxies, they could be compatible 
with an M82-like SED for z > 1.2 and 1^ > 4 10^^ Lq or with an Arp220-like SED (though 
fainter), if only their redshift is larger than z = 0.3. In any case, the FIRBACK galaxies are 
most probably LIRGs or ULIRGs and their large number density implied by the number counts 
implies a strong evolution of this population of galaxies with respect to IRAS, consistently with 
ISOCAM, although the 175 /im galaxies are most probably brighter than the 15 /xm ones. 

6 Conclusion 

The ISO extragalactic surveys have demonstrated the fundamental role of IR in the understand- 
ing of galaxy evolution. However, only ground-based follow-ups will allow us to understand the 
nature of the population of galaxies responsible for the strong evolution exhibited by the IR 
number counts. Are they particularly rich in metals ? Do they harbor active galactic nuclei ? 
We have a large enough sample of galaxies to be able to quantify the evolution of the luminosity 
function in the IR as a function of redshift. However, we are still far from understanding the 
origin of the huge cosmic IR background revealed by COBE, but SIRTF and FIRST are still 
to come... 
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